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ABSTRACT 


During  the  first  year  of  research  under  the  above  contract  the 
problem  of  the  stability  of  crystal  structures  for  the  heavy  rare  gas 
atoms  has  been  analyzed.  It  is  known  from  good  experimental  evidence 
that  neon,  argon,  krypton  and  xenon  crystallize  in  the  face-centered 
cubic  structure  (fee),  whereas  helium  under  pressure  crystallizes  in 
the  hexagonal  close-packed  form  (hep).  According  to  calculations  based 
on  twobody  forces,  however,  the  hexagonal  lattice  should  be  somewhat 
more  stable  for  all  solids  involved. 

The  underlying  hypothesis  for  the  analysis  has  been  that  the  cubic 
structure  becomes  more  stable  if  deviations  from  the  assumption  of 
twobody  forces  are  taken  into  account.  Since,  in  principle,  there  occur 
many  possible  types  of  such  manybody  interactions,  a  systematic  analy¬ 
sis  of  these  various  effects  was  first  undertaken,  based  on  a  cluster  ex¬ 
pansion  and  making  use  of  a  perturbation  treatment. 

Simultaneously,  a  number  of  pilot  calculations  was  carried  out, 
with  the  aim  of  acquiring  better  insight  in  the  different  characteristics 
of  the  various  possible  manybody  effects.  From  the  analysis  of  spectra 
from  trapped  free  radicals  in  molecular  crystals  it  had  been  found  that 
an  argon  matrix  behaves  to  some  extent  like  a  nitrogen  matrix.  To  pro¬ 
vide  a  basis  for  the  explanation  of  this  curious  phenomenon,  it  was  shown 
that  the  crystal  field  in  solid  argon  is  not  simply  the  sum  of  fields  from 
separate  atoms,  but  that  deviations  occur  in  the  form  of  higher  multipole 
moments  induced  by  electron  exchange  between  neighboring  atoms.  Such 
exchange  effects  imply,  for  example,  that  there  is  a  quadrupole  contribu¬ 
tion  arising  from  an  argon-diatomic  "molecule",  to  the  crystal  field.  It 
was  found  that  such  a  quadrupole  moment  is  not  small  compared  to  that 
of  a  nitrogen  molecule  and  that  it  may,  at  least  in  part,  account  for  the 
splitting  of  spectral  lines  emitted  by  excited  free  radicals  (nitrogen  atoms) 
in  an  argon  matrix. 


Also,  preliminary  calculations  were  performed  with  various  forms 
for  the  threebody  force  (different  angular  and  radial  functions).  The  dif¬ 
ferences  between  crystal  sums  for  the  hep-  and  fcc-lattices  increase  in 
general  as  the  forces  become  of  shorter  range,  and  increase  as  well 
with  higher  angular  dependence  (increasing  powers  of  cos  0,  where  0  is  the  angle 
between  the  lines  connecting  any  two  pairs  of  atoms  on  a  triangle).  The 
relative  difference  between  the  crystal  sums  may  be  of  either  sign  and 
is,  for  reasonable  analytic  threebody  potentials,  of  the  order  of  magni¬ 
tude  of  ten  to  twenty  percent. 

After  this  stage  of  the  research  had  been  completed,  the  general 
treatment  of  threebody  forces  between  rare  gas  atoms  in  first  order  of 
perturbation  theory  was  initiated.  To  avoid  difficulties  associated  with 
the  use  of  a  multipole  expansion  for  the  electrostatic  interactions  between 
the  atoms  (such  series  diverge  generally  for  all  interatomic  distances), 
the  integrals  were  evaluated  exactly  on  the  basis  of  a  model  calculation, 
in  which  the  electron  charge  distribution  of  the  atoms  is  approximated 
by  a  Gaussian  function  with  an  adjustable  parameter.  Values  for  this  pa¬ 
rameter  can  be  determined  empirically,  both  at  large  and  at  small  dis¬ 
tances  between  the  atoms. 

The  analysis  was  restricted  to  triangles  formed  by  a  central  atom 
and  its  twelve  nearest  neighbors  in  either  lattice  (fee  and  hep).  It  was 
found  that  the  relative  deviation  from  the  additive  first-order  energy  is 
negative  for  an  equilateral  triangle,  amounting  to  about  twenty  percent. 

For  a  120°-configuration  this  relative  difference  is  much  smaller  and 
always  positive,  whereas  for  a  linear  symmetric  array  the  effect  is  again 
large,  positive,  and  decreasing  practically  linearly  with  interatomic  dis¬ 
tance. 

A  comparison,  finally,  between  the  66  triangles,  formed  by  a  central 
atom  and  its  twelve  nearest  neighbors  in  the  hep-  and  fee -structures  (57 
are  the  same  between  the  lattices,  whereas  9  are  different)  revealed  that, 
although  this  threebody  effect  contributes  significantly  to  the  first-order 
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energy  of  the  crystals,  the  difference  between  the  two  structures  is  still 
too  small  to  account  for  the  observed  stability  of  the  face-centered  cubic 
structure. 

Therefore,  a  second-order  calculation  was  started  on  the  same  basis, 
leading  to  many  body  contributions  to  the  van  der  Waal  s -forces.  This  analy¬ 
sis  is  at  present  under  way. 


OUTLINE  OF  RESEARCH  PROJECT 


Concerning  properties  of  materials  in  general,  it  is  nowadays  re¬ 
cognized  that  no  essential  progress  is  possible  without  extensive  funda¬ 
mental  research.  In  the  past,  most  of  the  theoretical  investigations  in 
the  field  of  properties  of  solid  materials  have  been  directed  towards  an 
explanation  of  phenomena  associated  with  metallic  substances,  of  semi¬ 
conductors,  and,  to  a  lesser  extent,  of  ionic  crystals.  Compared  with 
this  group  of  materials,  the  theory  of  molecular  solids  (insulators)  has 
received  but  little  attention,  and  even  the  most  striking  facts  exhibited 
by  such  materials  often  lack  a  satisfactory  explanation.  It  has  become 
more  and  more  urgent  to  try  to  fill  this  gap,  since  molecular  solids  have 
become  of  increasing  importance  also  from  a  practical  point  of  view.  As 
an  example  it  might  be  mentioned  that  the  stability  of  the  close-packed 
cubic  crystal  structure  of  the  heavy  rare  gases  (neon,  argon,  krypton 
and  xenon)  cannot  be  explained  on  the  basis  of  existing  theories.  One  finds 
that  the  hexagonal  close-packed  crystal  form  (the  helium  structure)  should 
be  more  stable  than  the  cubic  structure,  and  this  result  is  remarkably 
insensitive  with  respect  to  the  precise  form  of  intermolecular  potentials. 
Concerning  this  problem  the  present  investigator  has  suggested  that  an 
explanation  may  be  found  on  the  basis  of  the  effect  of  manybody  forces  of 
the  van  der  Waals-type  (second  order  of  perturbation  theory)  on  the  stabi- 
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lity.  Preliminary  calculations  have  shown  that  this  explanation  may  be 
the  correct  one,  and  that  present  theories  of  molecular  solids  need  con¬ 
siderable  modification. 

The  crystals  of  nitrogen  and  carbon  monoxide  are  also  cubic 
close-packed  at  the  lowest  temperatures,  but  show  a  transition  to  the 
hexagonal  closest  packing  at  35°  K  and  61°  K,  respectively.  The  pre¬ 
sent  investigator  has  suggested  that  the  stability  of  these  crystals  and 
the  transition  mechanism  may  involve  orientational  interactions  (e.  g. 
those  due  to  permanent  electric  quadrupole  moments)  in  the  solids. 

The  cubic  structure  is  favored  by  quadrupole  interactions  in  the  solids, 
whereas  the  hexagonal  structure  is  preferred  by  anisotropic  components 
of  the  van  der  Waals  forces.  As  the  temperature  increases  and  herewith 
the  rotational  amplitudes  of  the  molecular  axes,  the  quadrupole  forces, 
being  of  first  order,  die  out  much  more  quickly  that  the  (second-order) 
van  der  Waals  forces.  The  results  obtained  so  far  suggest  strongly  that 
this  explanation  may  be  the  correct  one,  although  a  quantitative  theory 
is  still  lacking. 

Molecular  solids  have  recently  become  of  importance  also  as  ma¬ 
terials  in  which  free  radicals  can  be  trapped  and  stored.  The  results  ob¬ 
tained  in  this  field  provide  an  invaluable  experimental  background  for 
the  development  of  a  crystal  field  theory  of  molecular  solids.  This  in¬ 
vestigator  has  been  associated  with  the  Free  Radicals  Research  Group 
at  the  National  Bureau  of  Standards  in  Washington,  D.  C.,  and  has  per¬ 
formed  some  preliminary  calculations  which  seem  to  indicate  that  pre¬ 
sent  theories  of  molecular  solids  are  in  need  of  substantial  improvement 
in  order  to  yield  satisfactory  agreement  with  experimental  data. 
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SUMMARY  OF  RESULTS  OBTAINED  DURING  FIRST  YEAR 


The  results  obtained  so  far  may  be  divided  into  four  different  ca¬ 
tegories  : 


I.  Development  of  a  general  theoretical  method  to  classify 
all  possible  types  of  manybody  interactions  in  molecular  solids.  This 
analysis  is  based  on  the  method  of  Localized  Orbitals  in  Crystals.  It 
leads  to  a  double  series  expansion  for  arbitrary  properties  of  molecu¬ 
lar  crystals  :  one,  in  terms  of  linked  exchange-clusters  of  increasing 
size,  and  one,  in  terms  of  increasing  orders  of  perturbation  theory. 
(Special  Technical  and  Scientific  Reports  No.  1  and  2,  January  30  and 
July  31,  1961.) 

II.  Evaluation  of  the  triatomic,  first-order  manybody  effect. 

The  calculation  makes  use  of  a  Gaussian  effective-electron  model  for 
rare  gas  atoms  (Special  Technical  and  Scientific  Report  No.  2,  July  3  lf 
1961). 

III.  Application  to  the  stability  of  crystal  structures  of  the 
heavy  rare  gas  atoms,  for  a  central  atom  and  its  twelve  nearest  neigh¬ 
bors.  The  hexagonal  close-packed  and  face-centered  cubic  lattices  are 
compared. 

IV.  Application  to  the  crystal  field  in  molecular  solids,  in  con¬ 
nection  with  a  future  analysis  of  spectra  emitted  by  free  radicals  trap¬ 
ped  in  solid  argon.  The  exchange  quadrupole  moment  for  two  argon  a- 
toms  in  solid  argon  is  calculated  (Special  Technical  and  Scientific  Re¬ 
port  No.  2,  July  31,  1961). 


The  analysis  of  threebody  interactions  in  second  order  of  perturba¬ 
tion  theory,  again  based  on  a  Gaussian  model,  is  at  present  under  way. 

A  Special  Technical  and  Scientific  Report  No  4  (August  15,  1961)  was  is¬ 
sued,  containing  details  and  further  technical  information  concerning  the 
calculations  for  I  and  II. 
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ANALYSIS  OF  RESEARCH  PERFORMED 


A.  GENERAL 


The  following  sections  contain  a  detailed  discussion  of  results  ob¬ 
tained  on  the  problem  of  stability  of  crystals  of  the  heavy  rare  gas  atoms. 
Excluded  from  the  considerations  is  solid  helium,  since  this  element 
crystallizes  at  moderate  pressures  in  a  hexagonal  close-packed  lattice, 
which  result  is  predicted  even  by  a  simple  analysis  based  on  pair-inter¬ 
actions  between  the  atoms  in  the  crystal.  This  result  is  not  changed  by 
any  manybody  effect  of  the  type  to  be  discussed  below,  since  such  effects 
depend  essentially  on  the  amount  of  overlap  between  neighboring  atoms 
in  the  solid.  Due  to  the  small  size  of  the  helium  atom,  and  its  large  zero- 
point  energy  which  results  in  a  relatively  large  nearest  neighbor  distance, 
overlap  effects  are  entirely  negligible  in  this  case. 

The  analysis  of  manybody  interactions  in  molecular  solids  will  be 
based  on  the  method  of  Localized  Orbitals  in  Crystals  (also  called  the 
Crystal  Valence  Bond  method,  the  crystal  analogon  of  the  He itler- London 
method  for  molecules).  We  could,  of  course,  just  as  well  have  used  a 
Crystal  Orbital  method  (the  crystal  analogon  of  the  Molecular  Orbital  me¬ 
thod  for  molecules,  also  called  the  LCAO-,  or  Bloch-scheme,  or  the 
Tight-Binding  Approximation)  which  gives  results  identical  with  those  of 
the  Localized  Orbitals  in  the  case  of  closed  shells  of  electrons,  as  with 
the  rare  gas  atoms. 

Further,  we  are  considering  crystals  which  are  held  together  by 
relatively  weak  forces  of  the  van  der  Waals-type.  These  interactions 
are  very  much  weaker  than  those  involved  in  chemical  binding,  or  in 
ionic  crystals,  or  even  compared  with  the  cohesive  energies  of  metals. 

In  good  approximation,  the  physical  properties  of  these  systems  may  then 
be  described  by  means  of  perturbation  methods. 
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It  will  be  shown  below  that  any  property  of  the  crystal  (in  particu¬ 
lar,  the  static  lattice  energy)  can  be  expressed  in  terms  of  a  double  se¬ 
ries  expansion  :  one,  in  terms  of  contributions  due  to  linked  exchange- 
clusters  of  increasing  size,  and  one,  in  increasing  orders  of  perturba¬ 
tion  theory.  Molecular  crystals  may  thus  be  defined  by  the  requirement 
that  this  double  series  expansion  converges  rapidly  in  both  respects. 


B.  SYSTEMATIC  ANALYSIS  OF  MANYBODY  INTERACTIONS 
IN  MOLECULAR  SOLIDS 


An  analysis  is  undertaken  of  the  different  possible  types  of  simul¬ 
taneous  interactions  between  more  than  two  atoms  or  molecules  in  so- 
called  molecular  solids.  The  analysis  is  carried  out  on  the  basis  of  a 
double  series  expansion  :  1)  in  terms  of  linked  exchange-clusters  of 
increasing  numbers  of  atoms;  2)  as  a  series  in  increasing  orders  of 
perturbation  theory.  The  use  of  a  multipole  series  for  the  electrosta¬ 
tic  interactions  between  different  atoms  is  avoided  by  retaining  this 
interaction  in  unexpanded  form.  Instead,  an  effective-electron  model 
is  used  with  a  Gaussian  form  for  the  charge  distributions.  The  method 
is  illustrated  by  computing  the  exchange  quadrupole  moment  of  two 
argon  atoms  as  a  function  of  their  distance.  Calculations  by  P.  Rosen 
and  by  Shostak  for  first-order  interactions  between  three  helium  atoms 
are  extended  to  atoms  of  the  heavy  rare  gases.  It  is  found  that  the  re¬ 
lative  magnitude  of  this  manybody  effect  may  amount  to  twenty  percent 
of  the  first-order  interaction  energy.  Possible  implications  with  res¬ 
pect  to  stability  of  the  cubic  structures  of  heavy  rare  gas  crystals  are 
briefly  discussed. 
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INTRODUCTION 


The  possible  importance  of  simultaneous  interactions  between  more 
than  two  atoms,  molecules  or  nucleons  for  the  interpretation  of  proper¬ 
ties  of  compressed  gases,  liquids,  solids  and  nuclear  matter  has  from 
time  to  time  aroused  interest  in  the  literature.  In  molecular  physics, 
this  interest  arose  principally  from  a  possible  r61e  of  manybody  inter¬ 
actions  in  a  solution  to  the  problem  concerning  the  stability  of  the  ob¬ 
served  cubic  crystal  structures  of  the  heavy  rare  gases.  Calculations 
based  on  pair-interactions  had  revealed  that  a  hexagonal  structure 
should  be  somewhat  more  favorable,  in  contradiction  with  experiment. 

It  seemed,  therefore,  that  the  pair- assumption  (additivity)  of  interac¬ 
tions  between  rare  gas  atoms  could  not  be  reconciled  with  their  crystal 
structure. 

In  contrast  with  nuclear  problems  we  possess  for  molecular  sys¬ 
tems  complete  knowledge  regarding  the  origin  of  the  forces  between  the 
particles,  so  that,  in  principle,  it  should  be  possible  to  assess  the  si¬ 
gnificance  of  manybody  interactions  for  the  properties  of  molecular  so¬ 
lids,  for  example.  Yet,  the  explicit  calculations  which  have  so  far  been 
carried  out  were  either  based  on  approximations  which  clearly  ceased 
to  be  satisfactory  at  relatively  small  interatomic  distances,  or  they 
concerned  systems  of  little  interest. 

As  an  example  of  the  first  category,  we  mention  the  so-called 

triple-dipole  effect  for  three  rare  gas  atoms,  evaluated  by  Axilrod 
1  2 

and  Teller  ,  and  Axilrod  .  This  calculation  concerns  the  induced- 
dipole  interaction  in  third  order  of  perturbation  theory  between  three 
non-overlapping  distributions  of  charge.  As  such,  it  constitutes  simply 

1. B.M.  Axilrod  and  E.  Teller,  J.  Chem.  Phys.  _n»  299  (1943) 

2 .  B.  M.  Axilrod,  J.  Chem.  Phys.  17,  1349  (1949);  ^9,  719,  724  (1951) 
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an  extension  of  van  der  Waals-interactions,  as  calculated  by  London, 
to  third  order.  But  manybody  interactions  occur  already  in  lower  or¬ 
ders  of  perturbation  theory  (induced  by  electron  exchange);  in  addition, 
the  use  of  a  multipole  expansion  for  the  electrostatic  interactions  be¬ 
comes  increasingly  less  accurate  as  the  distances  between  the  atoms 

3  4  5 

decrease.  In  fact,  such  series  usually  do  not  converge  at  all  ’  ' 

6 

Tredgold  and  Ayres  have  modified  this  calculation  by  using  a  Gaussian 
distribution  function  for  the  negative  charge  cloud  of  each  atom,  and 
treating  the  overlap  of  charge  distributions  of  neighbouring  atoms  in  a 
classical  way.  This  adds  still  another  source  of  error,  since  at  small 
distances  the  interactions  are  mainly  of  exchange  type. 

The  same  type  of  third-order  effect  was  evaluated  independently 
7 

by  Muto  using  an  oscillator  model  for  the  atoms,  and  later  extended  by 

8  9 

Midzuno  and  Kihara  .  Bade  computed  the  dipole -dipole  part  of  the 
London- van  der  Waals  interaction  energy  in  the  general  order  of  per¬ 
turbation  theory,  without  exchange,  using  a  model  which  represents 
each  molecule  as  an  isotropic  harmonic  oscillator.  His  results  were 
applied  to  a  linear  chain  of  atoms  by  Bade  and  Kirkwood 3 4 5 6 7 8 9  10 11. 

On  the  other  hand,  first-order  forces  (exchange,  chemical  or 
valence  forces)  are  predominant  when  the  interatomic  distances  are 
small.  Margenau  had  already  pointed  out  that  they  are  of  a  manybody 


3.  F.  C.  Brooks,  Phys.  Rev.  86,  92  (1952) 

4.  A.  Dalgarno  and  J.  T.  Lewis,  Proc.  Roy.  Soc.  (London)  A69,  57  (1956) 

5.  L.  Jansen,  Physic  a  23,  599  (1957);  Phys.  Rev.  110,  661  (1958) 

6.  R.U.  Ayres  and  R.  H.  Tredgold,  Phys.  Rev.  100,  1257  (1955); 

Bull.  Ams.  Phys.  Soc.  L  Ser.  II,  292  (1956) 

7.  Y.  Muto,  Proc.  Phys. -Math.  Soc.  Japan,  62 9  (1943) 

8.  Y.  Midzuno  and  T.  Kihara,  J.  Phys.  Soc.  Japan  H,  1045  (1956) 

9.  W.  L.  Bade,  J.  Chem.  Phys.  27,  1280  (1957) 

10.  W.  L.  Bade  and  J.G.  Kirkwood,  J.  Chem.  Phys.  27,  1084  (1957) 

11.  H.  Margenau,  Rev.  Mod.  Phys.  J_l»  1  (1939). 
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type.  Such  an  effect  was  evaluated  for  an  equilateral  triangle,  and  a 

12 

linear  symmetric  array,  of  three  helium  atoms  by  P.  Rosen  ,  Shostak 
repeated  the  calculations  for  the  linear  array,  using  molecular  orbitals. 
Unfortunately,  helium  is  of  no  direct  interest  as  far  as  its  crystal  struc¬ 
ture  is  concerned,  since  this  element  crystallizes  in  a  hexagonal  latti¬ 
ce,  in  agreement  with  the  result  of  calculations  based  on  the  assump¬ 
tion  of  pair-interactions. 

It  is,  nonetheless,  of  interest  to  note  the  sign  of  the  relative 
threebody  interaction  energy  :  both  in  the  third-order  (Axilrod-Teller) 
and  the  first-order  (Rosen-Shostak)  effects  this  sign  is  minus  for  an 
equilateral  triangle,  and  plus  for  a  linear  array  of  atoms.  This  means 
that  for  an  equilateral  triangle  the  attraction,  respectively  repulsion  is 
weakened  compared  to  a  simple  sum  over  three  isolated  pairs,  whereas 
the  interactions  for  the  linear  array  are  stronger  than  the  pair-assump¬ 
tion  indicates. 

If  exchange  effects  are  taken  into  account,  then  manybody  interac¬ 
tions  occur  also  in  second  order  of  perturbation  theory,  i.  e. ,  in  the  sa- 

14  15 

me  order  as  the  van  der  Waals  forces  '  .  In  this  case  we  have  to  do 

with  a  triplet  of  atoms,  two  of  which  overlap,  whereas  exchange  with  the 
third  atom  is  not  considered.  Calculations  were  carried  out  on  the  basis 
of  a  Gaussian  model  for  the  electron  charge  distributions  of  the  atoms, 
for  the  dipole-dipole  and  dipole -quadrupole  components  of  the  interac¬ 
tions.  Again,  the  use  of  a  multipole  series  renders  the  results  unreliable 
for  small  interatomic  separations,  so  that  quantitative  conclusions  re¬ 
garding  crystal  stability  cannot  be  drawn. 


12.  P.  Rosen,  J.  Chem.  Phys.  21,  1007  (1953) 

13.  A.  Shostak,  J.  Chem.  Phys.  23,  1808  (1955);  see  also  P.O.  LOwdin, 
J.  Chem.  Phys.  19,  1570,1579^(1951)  (ionic  crystal) 

14.  R.  T.  Me  Ginnies  and  L.  Jansen,  Phys.  Rev.  101,  1301  (1956) 

15.  L.  Jansen  and  R.  T.  Me  Ginnies,  Phys.  Rev.  104,  961  (1956); 

R.  T.  Me  Ginnies,  Ph.  D.  Thesis,  Univ.  of  Maryland  (1957). 
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For  completeness,  we  mention  manybody  interactions  between 

atoms  or  molecules  which  are  originally  in  excited  electron  states. 

They  are  of  importance  at  very  high  temperatures  and  have  been  dis- 

1  fi 

cussed  by  Dahler  and  Hirschfelder  .  Such  effects,  together  with 

additional  complications  arising  in  case  the  charge  distributions  are 

15  16 

not  spherically  symmetric  *  ,  will  not  be  considered  in  this  paper. 

In  view  of  the  various  shortcomings  of  the  existing  calculations, 
it  seems  useful  to  devise  a  systematic  analysis  of  all  possible  types  of 
manybody  interactions,  and  to  carry  out  their  evaluation  in  a  consis¬ 
tent  manner,  avoiding  assumptions  which  invalidate  their  application  to, 
especially,  the  stability  of  rare  gas  crystal  structures.  Since  the  inter¬ 
actions  in  molecular  crystals  are  very  weak  compared  to  those  in  metals 
or  in  ionic  crystals,  we  will  use  a  perturbation  expansion  starting  from 
free-atom  wave  functions,  but  avoid  using  a  multipole  series  for  the 
perturbation. 


LINKED  EXCHANGE-CLUSTER  EXPANSION 


Starting  from  ground-state  free-atom  wave  functions  <Pa*<pb*  <pc»  • 
where _a,  b,  c, . . .  denote  atoms,  which  are  antisymmetric  with  respect 
to  exchange  of  electrons  of  the  same  atom,  we  write  the  zeroth  order 
wave  functions,  except  for  normalization,  of  the  assembly  of  atoms  as 
(Slater  determinant) 

<p(o)  =  2  <-l)X  £  <pa  <pc .  Cl) 

\ 


16.  J.S.  Dahler  and  J.  O.  Hirschfelder,  J.  Chem.  Phys.  25,  986  (1956) 
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an  antisymmetric  sum  of  products  of  atomic  wave  functions;  exchanges 
electrons  between  different  atoms,  or  it  denotes  a  combination  of  inter-  and 
intra-atomic  exchange. 

Relative  to  a  fixed  association  of  electrons  with  the  different  atoms, 
each  in  (1)  can  be  written  as  a  product  of  P's  involving  exchange- 
clusters  of  atoms.  For  example,  if  P^  denotes  exchange  of  electrons  1 
and  2  between  atoms  £  and  b,  2  and  3  between  atoms  Is  and  c,  5  and  6  bet¬ 
ween  in  and  n,  then  we  call  the  triplet  (ab  c)  a  linked  exchange -cluster; 
the  pair  (mn)  is  also  such  a  linked  cluster.  If  we  follow  a  perturbation 
method  for  the  evaluation  of  any  property  of  the  solid,  then  we  will  in 
principle  obtain  contributions  due  to  such  exchange-clusters  in  every  or¬ 
der  of  approximation.  These  contributions  can  be  ordered  according  to 
the  largest  cluster  which  they  contain  :  monatomic  terms  (without  exchange), 
diatomic  terms  (only  pair  exchange),  triatomic  contributions  involving  tri¬ 
plets  of  atoms  but  no  larger  clusters,  etc.  With  increasing  cluster-size  we 
describe  the  solid  in  terms  of  units  of  increasing  numbers  of  atoms  :  single 
atoms,  pairs,  triplets,  etc;  this  procedure  amounts  to  an  increasing  delo¬ 
calization  of  electrons  in  the  solid. 

In  this  way,  a  double  series  expansion  is  obtained  :  one,  in  terms 
of  exchange-clusters  of  increasing  size,  and  one,  in  increasing  orders 
of  perturbation  theory.  For  molecular  solids  this  procedure  may  be  ex¬ 
pected  to  yield  convergent  series  in  both  directions.  In  fact,  it  is  conve¬ 
nient  to  define  molecular  solids  by  the  requirement  that  both  the  cluster 
expansion  and  the  perturbation  series  converge  rapidly  for  aqy  property 
of  interest.  Such  a  solid  can  be  expected  to  show  molecular  characteristics 
to  a  large  extent. 

The  different  possible  types  of  manybody  interactions  can  now  be  in¬ 
dicated.  In  the  following  Table  I,  we  place  to  the  left  the  exchange -clusters 
of  increasing  size,  to  the  right  the  different  orders  of  perturbation  theory. 

A  straight  line  between  left-  and  right-columns  connects  a  cluster  term 
with  the  lowest  order  of  perturbation  theory  in  which  manybody  interactions 


occur. 
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Table  I.  Possible  types  of  manybody  interactions 


exchange -cluster 


order  perturbation  theory 


monatomic 

diatomic 

triatomic 


first 

second 

third 


If  we  use  two  sets  of  numbers  1,  2,  3, ... ,  to  characterize  a  type  of 

manybody  effect,  of  which  the  first  number  denotes  the  cluster-size,  the 

1  2 

second  the  order  of  perturbation  theory,  then  the  Axilrod-Teller  '  - 

14  15 

effect  is  of  type  1-3,  the  second-order  (van  der  Waals)  manybody  term  ' 

8  9 

of  type  2-2,  and  the  Rosen-Shostak  '  -calculation  of  type  3-1.  The  com¬ 
mon  twobody  interactions  may  also  be  indicated  in  this  manner  :  1-2,  i.  e. 
van  der  Waals  forces,  and  2-1,  i.e.,  repulsive  interactions  between 
closed  shells  at  small  interatomic  distances  (exchange  forces). 


A  GAUSSIAN  EFFECTIVE-ELECTRON  MODEL 


The  task  is  now  to  evaluate  these  different  types  of  marybody  inter¬ 
actions.  In  principle,  accurate  knowledge  of  atomic  wave  functions  is  re¬ 
quired  to  obtain  precise  results.  However,  if  we  have  in  mind  an  appli¬ 
cation  to  the  stability  problem  for  crystals  of  the  heavy  rare  gases,  then 
the  use  of  a  more  approximate  model  may  be  justified,  since  the  effect 
can  depend  sensitively  only  on  some  general  characteristics  of  the  wave 
functions  ("size"  of  the  atom,  plus  eventually  symmetry  properties). 

We  will  simplify  the  problem  by  means  of  the  following  assumptions  : 


.  15  - 


1)  only  single  interatomic  exchange  of  electron  pairs  is  taken 
into  account.  Effects  due  to  exchange  of  two  or  more  pairs 
of  electrons  between  the  same  two  atoms  are  therefore  not 
considered; 

2)  contributions  due  to  coupling  of  inter-  and  intra-atomic  ex¬ 
change  are  neglected; 

3)  the  electron  charge  distribution  of  an  atom  has  spherical 
symmetry. 

The  assumptions  1)  and  2)  must  be  applied  with  care.  In  many  cases 

the  overlap  integrals  between  the  atomic  orbitals  of  the  constituents  are 

not  small,  implying  that  higher-order  overlap  and  exchange  effects  may 
17 

not  be  neglected  .  In  the  case  of  rare  gas  crystals  effects  due  to  1), 

respectively  1)  and  2),  were  found  to  be  negligible,  except  for  very 

14  15 

small  interatomic  distances,  for  helium  and  neon  '  .  Validity  of 

assumption  3)  is  somewhat  in  doubt  for  atoms  of  the  heavy  rare  gases, 

18 

following  calculations  by  Linnett  and  Poe  ,  but  quantitative  results 
are  not  available. 

On  the  basis  of  the  above  three  assumptions,  the  following  model 
can  be  constructed.  We  sum  over  single-exchange  effects  between  all 
possible  pairs  of  electrons  (parallel  spins)  of  two,  or  three,  atoms.  The 
total  effect  is  then  replaced  by  exchange  between  one  "effective"  electron 
on  each  atom.  The  problem  becomes  thus  formally  the  same  as  that  con¬ 
cerning  hydrogen  atoms  with  parallel  spins.  The  charge  distribution  for 
the  effective  electron  is  chosen  to  be  of  Gaussian  form 


17.  cf.  P.  O.  LOwdin,  Adv.  in  Phys.  5,  1  (1956),  for  a  general  discus¬ 
sion  of  thk  difficulty. 

18.  J.W.  Linnett  and  A.  J.  Poe,  Trans.  Far.  Soc.  47,  1033  (1951). 

See  also  J.  Cuthbert  and  J.W.  Linnett,  Trans.  Far.  Soc.  54, 

617  (1958). 
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p  (r)  =  {p/n1^2  )3  exp  (-p2r2);  ...  (2) 

r_  is  the  distance  from  the  effective  electron  to  its  nucleus,  p  is  a  para¬ 
meter.  Values  of  pfor  atoms  of  the  heavy  rare  gases  are  determined 
empirically,  by  calculating  second-order  dipole  interactions  between  two 

Gaussian  atoms  at  large  distances  and  comparing  the  result  with  an  empi- 

15 

rical  potential  function  .  In  this  way  we  fit  the  outer  part  of  the  Gaussian 
distribution.  For  short-range  interactions  the  inner  part  of  the  charge  dis¬ 
tribution  also  plays  an  important  rdle,  so  that  it  becomes  necessary  to 
compare  such  interactions  between  Gaussian  atoms  with  the  repulsive  part 
of  an  empirical  potential  function.  Such  a  comparison  will  be  carried  out 
in  a  following  section. 


EXAMPLE  ;  THE  EXCHANGE  QUADRUPOLE  MOMENT 
OF  A  PAIR  OF  ARGON  ATOMS 


To  illustrate  the  use  of  the  Gaussian  model,  we  calculate  the  elec¬ 
tric  quadrupole  moment  of  a  pair  of  argon  atoms,  a  distance  R  apart. 

The  exchange-repulsion  between  closed  shells  will  give  rise  to  a  small 
net  positive  charge  in  the  region  between  two  argon  atoms,  and  a  negative 
charge  elsewhere.  This  effect  gives  rise  to  exchange  multipole  moments 
of  even  order,  of  which  we  determine  the  first  one,  an  electric  quadru¬ 
pole  moment. 

Consider  two  Gaussian  atoms  at  interatomic  distance  R.  The 
charge  distribution  of  this  pair  has  cylindrical  symmetry  about  an 
axis  z,  coinciding  with  the  direction  of  R.  For  such  cylindrically 
symmetric  distributions  of  charge  we  can  define  scalar  multipole 
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5 

moments  of  arbitrary  (positive)  order  ;  the  first  non-vanishing  mo¬ 
ment  is  in  this  case  a  quadrupole,  defined  by 

Q  =  <  2  e.z.2  >  -  <  2e.x.2>  ;  (3) 

i  i 

where  the  x-axis  is  perpendicular  to  z,  but  otherwise  arbitrary.  For 
two  Gaussian  atoms,  s^and  b,  we  have  two  electrons,  ^and  2,  and 
the  analogy  with  the  case  of  two  hydrogen  atoms  with  parallel  spins 
makes  it  possible  to  write  the  zero-order  wave  function  as 

,<0>  =  (  3(i.  J  )  •  <4> 

2 

where  <p  =  p  ,  the  Gaussian  charge  density  of  an  atom,  and  a  is 
the  overlap  integral,  a  -  f  p  (p^  dt  .  The  direct  terms  in  2 
give,  of  course,  no  contribution  to  Q  or  to  any  other  permanent 
multipole  moment,  so  that  we  are  left  with  the  exchange  terms  (omit¬ 
ting  normalization  )  -2  <p  a(l)  <^(1)  9a(2)  (^(2). 

We  denote  the  z-coordinates  of  the  position  vectors  of  electrons 
1  and  2,  with  respect  to  the  center  of  the  line  connecting  the  two  nu¬ 
clei,  by  Zj  and  z2,  respectively.  Then  we  obtain  : 

<2eiZl2>  -  ^  4e  1 3  A  -  4e  (R/2)2  A  2  |  /2  ( 1  -  A2), 

where 

h  sf  z2(Pa9b  dT; 


(5) 
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<  2  e.x.2  >  =  4  e  If  /2  (1  -A  2) 

with 

l2  =  /*2  a  *Pb  *  ‘  (6) 

The  expression  for  the  scalar  quadrupole  moment  of  the  pair 
(ab)  of  argon  atoms  becomes  with  (3),  (5)  and  (6), 

Q  =  <2  e  z  2  >  -  <  2  ejX  2  >  = 

i  1  1  i  1  1 

=  1 4  e  a(I3  -  I2  )  -  4  e  a2(R/2)2  J  /2  (1  -  a2)  (7) 

For  a  Gaussian  distribution  I  ^  =  i2,  so  that  then 
2 

q  =  _2e  - h—  (R/2)2.  (8) 

1  -  A 

This  quadrupole  moment  is  the  same  as  that  caused  by  an  effec¬ 
tive  negative  point-charge  -<6e)  =  -e  — - --  at  nucleus  a,  the  same 

1  -  AZ 

charge  at  nucleus  b,  and  a  positive  charge  2(  6  e)  at  the  center  bet¬ 
ween  the  two  nuclei. 
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2  2 

Values  of  A  /(I  -  A  )  can  be  determined  with  the  help  of  those  of 

the  parameter  J3  in  the  Gaussian  distribution  function,  for  the  heavy  rare 

15  2  2 

gases  .  For  nearest  neighbors  in  solid  argon  A  /(I  -  A  )  =  0.06,  and 

R/2  *  1.92  &  .  The  resulting  value  for  the  scalar  quadrupole  moment  of 

a  pair  of  nearest -neighbor  argon  atoms  is  then 

Q  =  2.  14  x  10”^®  e.s.u  =  0.45  x  10"  e  cm^, 

i.  e. ,  a  value  of  the  same  order  of  magnitude  as  for  a  nitrogen  molecule. 
This  large  value  is  not  due  to  a  large  displacement  of  charge,  but  in  the 
first  place  to  the  large  value  of  the  distance  between  two  argon  atoms, 
which  enters  as  the  square  in  the  expression  for  the  quadrupole  moment 
(compare  a  nitrogen  molecule  :  R/2  =  0.55  $1 ). 

In  view  of  the  approximations  which  are  inherent  in  this  type  of 

model,  quantitative  conclusions  do  not  seem  to  be  justified.  The  result 

indicates,  however,  that  effects  due  to  diatomic  exchange-clusters  of 

heavy  rare  gas  atoms  may  be  significant  for  the  determination  of  some 

* 

solid  state  properties  . 


THE  TRIA  TOMIC  -  FIR  ST  ORDER  (3-1)  MANYBODY  EFFECT 


The  task  is  now  to  evaluate  the  different  possible  types  of  many- 
body  interactions  in  first,  second  and  third  orders  of  perturbation  theory. 
The  remainder  of  this  paper  is  devoted  to  a  calculation  of  the  3-1-effect, 


This  quadrupole  effect  may  possibly  be  observed  by  means  of  spectra 
emitted  by  trapped  radicals  in  rare  gas  matrices,  e.g.  from  excited 
nitrogen  atoms  or  excited  alkali  atoms  in  solid  argon. 
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involving  a  cluster  of  three  atoms,  in  first  order  of  perturbation  theory. 

It  was  mentioned  earlier  that  P.  Rosen  and  Shostak  computed  such  inter¬ 
actions  for  three  helium  atoms.  Helium  is,  however,  of  no  direct  interest 
as  far  as  its  crystal  structure  is  concerned,  and  a  direct  extension  of  the 
methods  used  by  these  authors  to  the  heavy  rare  gases  would  be  very  com¬ 
plicated,  since  detailed  knowledge  of  the  wave  functions  is  required.  We 
will,  therefore,  make  use  of  the  Gaussian  model  ,  one  effective  electron 
per  atom,  with  a  charge  density  given  by  (2),  and  variable  0. 

Consider  such  a  triangle  (a  be)  of  atoms.  The  zero-order  wave 
function  is  (Slater  determinant) 

¥(o)  «  - - 1 - —  Det  j9a(l)<Pb(2)<Pc(3)|  ,  (9) 

'l  3  !  (1-  A2abc) 


with 


2  2  2  2 
A  abc  =  A  ab  +  A  ac  +  A  be  -  2  A  ab  Aac  Abe, 


(10) 


The  perturbation  Hamiltonian,  H'abc  ,  can  be  written  as 


H'  .  *  H'  +H*  +  H'  . 

abc  ab  ac  be 


For  the  first-order  energy  we  have  to  evaluate 


E1  •  <H'abc>  *  <H'ab  >  +  <H’.c>  +  <H'bc> 


7 
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After  substituting  the  expressions  for  H'ab  and  into 

the  following  result  is  obtained  : 


2  1  1  -  J  <Aac  +  A\c  > 

<H'ab  >/e2  =  — - 2  - - - - —  G__^  + 


R 


ab 


1  -  A 


abc 


aa(b) 


,  „  Aab  "  Aac  Abc  ^  ^  Aac  “  Aab  Abc  ^ 

+  2  - g -  Gab(a)  +  - ; - T -  Gac(b)  + 


1  -  A 


abc 


1  -  A 


abc 


+  Abc  ~  Aab  Aac  1  /A  A  )  + 

2  Gbc(a)  “  ?  (Aabab  "  Aaabb) 


1  -  A 


abc 


1  -  A 


abc 


ac 


+  <- - 2 

1  -  A 


>  <Aabbc  -  Aabcb>  +  <  :  *V  >  (A..bc  -  «» 


abc 


1  -  A 


abc 


where  the  symbols  G  and  A  are  abbreviations  for  the  following  integrals 


Gab(c) 


J  *a  fb 


:  Gaa(b)  = 


dx 


etc. , 


f  r  =  distance  between  an  electron  and  nucleus  of  atom  c, 
c 

)  rfa  =  distance  between  an  electron  and  nucleus  of  atom  1)  . 
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abac 


if 


<pa(i)  <pb(2)  9a(l)  <pc(2) 


dx. 


; 


12 


aabc 


if 


<PaU>  <pa(2)  cpb(l)  <pc( 2) 


dTl  *2  5 


12 


etc.,  rl2  =  distance  between  electrons  1  and  2. 


The  first-order  energy  between  a  and  b,  without  c  present,  is 


<  H'ab>  '  /e 


(o),  2  „  1 


Rab  1  “  *  ab 


5  t  2Gaa(b)  *  2AabGab(a) 


'  (Aabab  "  Aaabb)  ^  * 


(12) 


2 

The  threebody  component  of  <H'ab  >  /e  ,  namely 
{<  H'ab  >  -  <H'ab><0’}/e2. 

is  obtained  by  subtracting  (12)  from  (11).  We  do  not  need  any  further 
equations,  since  <H'  >,  for  example,  is  obtained  from  ( 1 1)  by  re- 

flv 

placing  b^  by  £  . 

The  integrals  G  and  A  have,  for  Gaussian  distributions  of  charge, 
19 

been  evaluated  by  Boys  .  For  the  atomic  wave  functions  we  take 


19.  S. F.  Boys,  Proc.  Roy.  Soc.  A200,  542  (1950). 
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9  =  =  <0/Il1/2)3/2  exp  (-  ^r2/2) 


It  is  convenient  to  use  the  following  abbreviation 


vr 


F(x)  = 


l=~  f 


dy  ; 


also,  we  denote  by  the  length  of  the  line  connecting  the  nucleus 

of  atom  a_  with  the  center  of  the  line  ,  etc.  The  following  table 
contains  a  list  of  the  integrals  occurring  in  the  equation  for  E^. 


Table  II.  List  of  integrals  for  E^ 


Integral 

A2ab 

A2ac 

A2bc 

A2abc  «A2ab  +  A2ac  +  A2bc  - 


Value  for  Gauss-distribution 
exp  (  -  f?R2  ab/2) 

exp  (  -02R2  ac/2) 

exp  (  -  02R2  bc/2) 

2Aab  Aac  Abe 
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G 

G 

G 

G 

A 

A 

A 

A 

A 

A 


aa(b)  =  Gbb(a) 

(2p/fF)  F(02R2ab) 

ab(a)  =  Gab(b) 

(2p/fif)^bF(  f32R2ab/4) 

bc(a) 

(2P/fn’)AbcF(p2R2a(bc)) 

ac(b) 

<2P/nry>  F(  (32R2b(ac)) 

abab 

P  \|T/  D  F  (  p2R2ab/2) 

aabb 

P  fi/fl  A2ab 

abac 

Pf2?n  Abc  F(p2R2a(bc)/2) 

abcb 

P(2?n  Aac  F(p2R2b(ac)/2) 

aabc 

Aab^cF(P2R2bc/8) 

abbc 

Pf*7n  AabVF(p2R2ac/8) 

MACHINE  CALCULATIONS 


The  computations  were  carried  out  on  the  digital  computer  Bull 
Gamma  AET  (word-length  :  12  digits;  memory  access  :  172  microseconds; 
type  of  storage  :  magnetic  drum  of  8192  memories).  After  calculating  the 
arguments  of  the  F-  and  a  -functions  for  various  triangular  configurations, 
these  functions  were  evaluated  by  interpolation,  using  "Tables  of  the  Error 
Function  and  its  Derivative11  (Nat.  Bureau  of  Standards,  Washington,  1953). 
The  interpolation  was  carried  to  the  fifth  term  of  a  Taylor  expansion;  this 
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appeared  sufficientfor  the  accuracy  desired.  The  intermediate  results 
were  transmitted  by  punchcards  to  the  main  program,  i.  e.,  the  compu¬ 
tation  of  <'H'a|3c>I  *  ^  their  difference,  and  the  relative 

threebody  component  of  the  first-order  interactions,  for  values  of  pR 
from  1.5  to  3.5,  in  intervals  of  0.  1,  for  each  triangular  configuration. 

For  values  of  pR  larger  than  about  3  the  results  are  somehwat 
irregular,  since  the  different  terms  in  equations  (11)  and  (12)  are  then 
very  small,  and  the  sums  and  differences  of  a  large  number  of  such 
terms  have  to  be  determined.  The  program  was  devised  in  such  a  way 
that  different  triangles  require  only  slight  modifications. 


RESULTS  FOR  SOME  SPECIAL  TRIANGLES 


The  threebody,  first-order  interaction  energies  were  evaluated  for 
a  number  of  special  triangles,  among  others,  an  equilateral  triangle,  a 
120°- symmetric  array  and  a  linear  symmetric  array.  These  three  types 
of  triangles  occur  among  those  formed  by  a  central  atom  and  two  of  its 
twelve  nearest  neighbors  in  the  cubic,  and  hexagonal,  close-packed  lat¬ 
tices.  The  relative  threebody  contribution  is 


<H’abc> 


-  <H'abc> 


<o) 


W 


<H'  K  > 
abc 


(13) 


where 


<  H'  .  >  =  <H'  .  +  <H'  +  <H'.  >*°* 

abc  ab  ac  be 


,  a  sum  of  first 
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order  energies  over  three  isolated  pairs  of  atoms.  The  relative  three- 

body  component  (13)  is  a  function  only  of  the  dimensionless  parameter 

(3  R,  where  R  is  the  smallest  distance  between  neighbors  in  the  triangle. 

Calculations  were  performed  for  the  three  types  of  triangles  with  PR  =  1.5 

1.  6,  . . . ;  3.  5,  in  intervals  of  0. 1.  This  range  covers  the  PR-values  for 

neon,  argon,  krypton  and  xenon,  which  are  3.44,  2.40,  2.  10  and  1.99, 

15 

respectively  .  The  results  are  collected  in  the  following  Figure  1. 


Fig.  1  :  Relative  threebody,  first-order  interactions  for  an  equilateral 
triangle,  a  120°-symmetric  and  a  linear  symmetric  array,  as  a  function 
of  PR. 


The  main  characteristics  of  the  threebody  interactions  are  : 

1)  Equilateral  triangle  .  The  relative  threebody  interaction  is  practically 
independent  of  P  R,  i.  e.  of  the  size  of  the  triangle  (at  constant  P  ). 

It  amounts  to  about  20  percent  and  its  sign  is  minus,  meaning  that  the 
total  first-order  interaction  is  weaker  (less  repulsive)  than  obtained 
on  the  basis  of  an  additive  sum-over-pairs.  The  sign  agrees  with  the 
result  obtained  by  P.  Rosen  for  three  helium-atoms; 

2)  120°-symmetric  array.  For  this  configuration  the  relative  threebody 
component  is  only  at  most  5  percent  of  the  additive  first-order  interac¬ 
tions  and  its  sign  is  plus  ,  signifying  slightly  enhanced  repulsion 
compared  with  the  sum  over  three  isolated  pairs; 
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3)  Linear  symmetric  array.  The  relative  threebody  interaction  is  very 
nearly  decreasing  linearly  with  increasing  pR,  from  20  percent  at 
PR  =  1.5  to  a  few  percent  at  pR  =  3.  Its  sign  is  plus,  which  means 
that  the  effect  enhances  the  repulsion  between  the  three  atoms,  com¬ 
pared  with  an  additive  sum-over-pairs.  The  sign  agrees  with  that  ob¬ 
tained  by  P.  Rosen  and  by  Shostak  for  three  helium  atoms. 


VALIDITY  OF  GAUSSIAN  MODEL  AT  SMALL  DISTANCES 


At  small  interatomic  distances  the  inner  part  of  the  charge  distri¬ 
bution  plays  an  important  rOle  in  the  determination  of  the  repulsive  in¬ 
teractions  between  closed  shells  of  electrons.  Since  the  values  for  the 
parameter  p  in  the  Gaussian  distribution  were  determined  empirically 
from  a  comparison  with  long  -  range  interactions,  it  is  of  interest  to 
check  whether  or  not  p  changes  appreciably  with  decreasing  distance 
between  the  atoms.  To  this  end  we  use  the  repulsive  part  of  a  modified 
Buckingham  potential. 


E(r)  «  - - -  (6/a)  exp  a  (1  -  R/R_), 

1  -  (6/a)  0 

with  values  of  e,  a  ,  and  R  for  neon,  argon,  krypton  and  xenon  as 

°  20 

given  by  Hirschfelder,  Curtiss  and  Bird  as  an  empirical  basis,  and 


20.  Hirschfelder, 
Liquids  (John 


Curtiss  and  Bird,  Molecular  Theory  of  Gases  and 
Wiley  and  Sons,  Inc.  New  Ifork,  1954)  Chap.  12. 
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compare  the  results  with  <H'  eq.  (12),  taking  for  0 the  values 

for  the  heavy  rare  gases,  and  varying  R.  The  results  of  this  compari¬ 
son  are  given  in  the  following  Figure  2. 


Fig.  2  :  Comparison  between  repulsive  part  of  modified  Buckingham 
potential  and  <H'ak>  ^  for  a  Gaussian  distribution,  for  the  heavy 
rare  gases,  as  function  of  interatomic  distance  R. 


It  is  seen  from  Figure  2  that  the  Gaussian  model  yields  excellent 
values  for  the  repulsive  potential  between  two  neon  atoms,  but  that 
with  increasing  atomic  weight  the  two  curves  deviate  more  and  more. 
This  implies  that  0  for  argon,  krypton  and  xenon  increases  some- 
what  as  R  decreases,  indicating  a  certain  "quenching"  of  the  three- 
body-effect.  It  should  be  remembered,  however,  that  with  decreasing 
R  multiple-exchange  effects  become  more  and  more  important, 
rendering  the  single-exchange  model  less  and  less  valid.  For  neon, 
argon,  krypton  and  xenon  the  nearest  neighbor  distances  in  the  crystals 
are  3.20  ft,  3.  84  ft  ,  3.94  ft  ,  and  4.37  ft,  respectively. 


CONCLUSIONS 


1)  It  is  possible  to  classify  the  various  manybody  interactions  bet¬ 
ween  atoms  in  "molecular  solids"  on  the  basis  of  a  double  series  ex¬ 
pansion;  one,  in  terms  of  linked  exchange-clusters  of  increasing  size, 
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and  one,  in  increasing  orders  of  perturbation  theory.  Molecular 
solids  may  then  be  defined  by  the  requirement  that  both  expansions 
converge  rapidly; 

2)  With  the  help  of  a  Gaussian  effective-electron  model,  the  triatomic- 
first-order  threebody  interactions  given  by  P.  Rosen  and  Shostak  for 
helium  atoms,  may  be  extended  to  the  heavy  rare  gases.  Both  for  an 
equilateral  triangle  and  a  linear  symmetric  array  the  sign  of  the  rela¬ 
tive  threebody  interaction  agrees  with  previous  results  for  helium; 

3)  For  a  single  triangular  configuration  the  relative  threebody  interac¬ 
tion  may  amount  to  20  percent  of  the  additive  first-order  energy; 

4)  Values  for  the  parameter  (3  of  the  Gaussian  distribution,  obtained 
from  a  comparison  with  long-range  interactions,  give  excellent  values 
for  short-range  forces  between  two  neon  atoms.  For  argon,  krypton 
and  xenon  0  appears  to  increase  somewhat  with  decreasing  interatomic 
distances. 


C.  THREEBODY  INTERACTIONS  BETWEEN  RARE  GAS  ATOMS 
IN  FIRST-ORDER  PERTURBATION  THEORY 


Compilation  of  Results 


General 

The  method  used  for  the  evaluation  of  threebody  interactions  in 
first-order  of  perturbation  theory  is  that  developed  in  section  B.  For 
completeness,  some  of  the  general  equations  will  be  repeated  from  this 
section. 


-  30  - 


The  triatomic -first -order  (3-1)  effect 

To  evaluate  the  3-1  effect  for  a  triangle  of  rare  gas  atoms,  use 
is  made  of  a  Gaussian  effective -electron  model.  The  total  effect  for  all 
electrons  of  the  three  atoms  is  replaced  by  that  due  to  one  "effective" 
electron  on  each  atom.  The  charge  distribution  of  the  effective  electron 
is  chosen  to  be  of  Gaussian  form  : 

P  (r)  =  (0/n  ^  )3  exp  (-  p2r2)  (14) 


where  0  is  a  parameter  which  can  be  determined  empirically. 

Consider  such  a  triangle  of  three  Gaussian  atoms  £ ,  t>  and  c, 
each  with  one  electron  (spins  parallel).  The  electrostatic  interaction 
between  the  three  atoms  can  be  written  as  : 


H*abc  "  H'ab  +  H'ac  +  H'bc  ' 


(15) 


whereas  the  first-order  energy  has  the  form  : 


<H,abc>  *  <Hab>  +  <H'ac>  +  <H'bc>  • 


(16) 


with  : 


<  H'abc  >  ■  ffj  *o*H'abc  *0  dV  dV  d*3-  etC' 

The  zero-order  wave  function  ¥  is  a  Slater  determinant  : 

o 


*0  =  ir{Tu~2~  ~)  D6t  fca(1W2><Pc<3>l  '  <17> 

’  *  ~  ^  abc' 

where  <p  ,  f .  and  ®  are  (antisymmetrized)  ground- state  wave  functions 
a  d  c 

for  the  atoms  a  ,  la  and  £ ,  and 

^  abc  '  ^  ab  +  ^  ac  +  ^  be  2^ab^ac^bc  ' 
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with 

Aab  =  J  tbe  over*aP  integral  between  a  and  b,  etc. 

If  we  denote  by  <H'abc  ^  ^  the  simple  sum  of  first-order  inter¬ 
actions  for  three  isolated  pairs  of  atoms,  then 

<H'abc  >  -  <  ^'abc  ^  ^  S^ves  the  threebody  component  of  the 
first-order  energy.  It  is  only  necessary  to  compute  <  H'  .  >  and 

/  \  3D 

<H'ab>  '  ,  since  the  remaining  terms  are  obtained  by  changing  in¬ 
dices.  The  results  are  : 


<H'ab> 


-  2 


R 


ab 


1  '  I  (AL  +  Abc  > 


1-  A 


abc 


Gaa(b)  + 


+  2  Ab-f  A.a9_*b£.  G 

1-Aabc 
^c  “  Aab  Aac 


A  ac  ~  Aab 


ab(a) 


1-  A 


1-  A 


Gbc(a)  + 


abc 

(A 


Gac(b)  + 


abc 


1- A 


abab  "  ^aabb^  + 


abc 


ac 


1-  A 


^c 

2“”  ^abbc  -  ^abcb^  +  "  ~2  ^aabc  "  "^abac^  * 


abc 


1-  A 


abc 


(18) 


where  the  symbols  G  and  A  are 
?a  ?b 


ab(c) 


*  /  r 


d  i  ; 


c 

r  =  distance  between  an 
c 

r,  =  distance  between  an 


abbreviations  for  the  integrals 
r  Va  ‘Pa 

Gaa(b>  ’  / -  *  ;  etc- 

J  rb 

electron  and  nucleus  of  atom  c, 
electron  and  nucleus  of  atom  b. 
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abac 


aabc 


■  If 
If 


<PaU><Pb(2>9a(1)<Pc(2) 

ri2 

<ftt(1)9a(2><Pb(1)  9c<2) 


d  Tj  d  t2  ; 


d  Tj  d  Tg  ;  etc. 


12 


r^2  =  distance  between  electrons  1  and  2. 

The  first-order  energy  between  a  and  b,  without  £  present,  is  : 


<  H'ab  >  “V 

Rab  1_  Aab 


^2<"*aa(b)  ”  2^ab^ab(a)~  ^abab“^aabb^ 

(19) 


Since  the  <p  's  are  real,  we  can  take  : 

<p  (r)  =  p  (r)1^  =  (  0/n  l/Z)  3/2  exp  <-p2  r2/2). 


(20) 


It  is  also  convenient  to  introduce  the  abbreviation  : 

,  ^  2 
F(z)  *  •— =—  f  e"  “  do  ; 

{  z  J 

o 

F  is  related  to  the  error  function  erf  (y)  or  H(y), 

y 


I 


.-a2 


d  a 


(21) 


as  follows  : 


F(y2)  -  —  f  e"  “2  da 

V  J 


2y 


H(y), 


(22) 


The  following  Table  II  contains  the  different  types  of  integrals, 
with  their  values  for  a  Gaussian  distribution  function. 
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TABLE  II 


Integral 


List  of  Integrals  for  the  3-1  Effect 

Value  for  Gauss  distribution 


A2h 

ab 

exp(-  p2R2ab/2) 

4 

exp(-  P2R2ac/2) 

4 

e*p(-  P2R2bc/2) 

4c  ' 

2  2 

A  ,  +A  + 
ab  ac 

4  -  2i.b 

A  A. 
ac  be 

Gaa(b) 

=  Gbb(a) 

(2  Pfl¥)  F(  ^R2ab) 

Gab(a) 

*  Gab(b) 

(?P$T>  F(p2R2ab/4) 

Gbc(a) 

(2P^bcF(p2R2a(bc)) 

Gac(b) 

<2 

Aabab 

P^T*  F(p2R2ab/2) 

Aaabb 

Pj2/^ 

Aabac 

P'T2^abcF(P2R2«(bc)/2» 

Aabcb 

Pf*A  AacF(P2R2b(ac)/2) 

Aaabc 

P^A  ^\cTi^R\c/S) 

Aabbc 

P|V*  ^b^FtpV^/e) 

Ra(bc) 

*  distance 

between  a 

and  the  middle  of  the  line  Rfac;  etc. 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


MACHINE  CALCULATIONS 


The  calculations  were  carried  out  on  the  digital  computer 
BULL  GAMMA  AET.  The  general  machine  features  are  : 

word-length  :  12  digits  (binary  commands) 

instructions  per  word  :  3 

operation  times  :  memory  access  172  p  s 

type  of  storage  :  magnetic  drum  of  8192  memories 

input -output  : 

punch  cards  equipment  :  speed  175  cards  per  minute 
cards-reader  equipment  :  speed  150  cards  per  minute 
printer  :  speed  150  lines  per  minute. 

programming  :  in  machine  language  or  in  symbolic 

programming  system  (AP2). 


Evaluation  of  eqs.  (18)  and  (19) 

The  threebody  first-order  effect  was  evaluated  for  several  diffe¬ 
rent  types  of  isosceles  triangles  (abc). 
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All  distances  are  made  dimensionless  through  multiplication  with 
the  Gaussian  parameter  p  .  We  need  the  following  quantities  : 


PRab  = 

8R 
^  ac 

e^c  = 

u^ab 

PRa(bc) 

*  PRab  ' 

u2 

1  ’  T~  =  vpRab 

{3Rab 

\|  n  2  „ 

pRb(ac) 

"  — 

U+2u  =wpRab 

where  u  =  2  sin  8/2,  v  =  lj  1  --37—  »  w  *  (  |  1  +  2  u2)/2. 

The  quantity  x  *  pR&b  can  be  taken  as  parameter,  in  terms  of 
which  : 

<  H'  .  >  /|3  e2  and  <H'  >  (o)  /Pe2 

«*d  gb 

are  given  for  any  particular  triangle.  A  range  of  values  for  x  was  chosen 
between  1.5  and  3.5,  in  intervals  of  0.  1. 

The  equations  to  be  evaluated  are  polynomials  consisting  of  expo¬ 
nential  terms  and  error  functions  : 


Exponentials  : 

(A) 

4b  =  exp  ("  x2/2> 

(B) 

4c =  exp  (_  x2/2) 

(C) 

=  exp  (-  ux2/2) 

^bc  "  ^ab  +  ^hc  +  ^bc  2  Aab  ^c^bc  * 
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Error  functions  :  The  G-  and  A-integrals  assume  the  following  simple 
forms  : 

(D)  Gaa(b)/P  s  (2A*>  F(x2)  =  °/x)  H(x) 

(E)  Gab(a)/p  =  (2/fi)AabF(x2/4)  =  (2/x)  H(x/2)^fa 
(F>  Gbc(a)^  =  <2 ^bcF<v2x2)  =  (1/vx)  H(vx)Abc 
<G>  Gac(b)/P  =  <^)AacF(w2x2)  =  (l/wx)  Htwx)^ 

(H)  Aabab/P  =  (1^T)  (2/^>  E(x2/2)  =  (l/x)  H(x/f2) 

(J>  Aaabb/P^  *ab 

(K)  Aabac/p  =  (1/^2>  (2/T^bcF(v2x2/2)  =  (l/vx)  HCvx/fT)^ 

(L)  Aabcb/p  *  (1/'T2)  (2/fi0^acF(w2x2/2)  *  (l/wx)  H(wx/f2)Aac 

(M)  Aaabc/P  =  (1/'T2>  (2/f^)AabAacF(u2x2/8)=(2/ux)H(ux/2f2)AabAac 

(N)  Aabbc/P  “  (1/'fT)  <2/f^)AabAbcF(x2/8)=(2/x)H(x/2'T2)AabAbc 

Symbolic  Programs  : 

It  appeared  essential  to  employ  a  floating  point  program.  The 
programming  was  carried  out  in  small  part  in  machine  language  and  the 
rest  in  symbolic  language  which  reduced  the  work  considerably.  The 
computations  can  be  subdivided  into  the  following  independent  parts  : 
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a)  Calculation  of  arguments  of  the  G-  and  A-integrals.  This  very 
short  program  was  primarily  intended  as  a  check  on  the  sym¬ 
bolic  program  AP2. 


b)  Calculation  of  the  exponentials  and  error  functions  : 

x 

H'(x)  =  (2/nn  exp(-x2);  H(x)  =  (2/jTn  J  e"a  da. 

o 

First,  H'(x)  was  computed  by  means  of  a  subprogram  (Bull)  with  an 
error  less  than  1  in  the  10th  decimal  place.  Then,  H(x)  was  evaluated  by 
direct  integration  of  H'(x)  with  the  help  of  Simpson's  method  : 


H<Vl>  =  H(Vl>  +  4  t  +  4  H'(xn>  4  J  • 

with  an  error  «  H(V)h5/90; 


*2m 


J  H'(x)dx  =  -g-  [  H'(o)  +  4H'(1)  +  2H'(2)  +  4H'(3)  +  . .  .  + 

+  4H  2m-1  +  H  2m  ^  ' 
with  an  error  <  mH^V^h^/90. 


This  method  had,  however,  to  be  abandoned,  since  the  necessary 
high  precision  could  not  be  reached. 


A  second  program  was  based  on  interpolation  between  values  given 
in  "Tables  of  the  Error  Function  and  its  Derivative11  (National  Bureau  of 
Standards,  Washington  D.  C. ,  1953).  These  tables  give  values  for  H'(x) 

_4 

and  H(x)  to  15  decimals  in  intervals  of  10  for  x  between  0  and  1,  and 
-3 

in  intervals  of  10  for  x  between  1  and  5.  6. 


It  appeared  more  advantageous  to  employ  a  high-order  interpola¬ 
tion  formula  and  fewer  tabulated  values  than  inversely.  This  avoided  the 
necessity  of  storing  long  lists  of  tables  in  the  machine.  A  Taylor-expansion 
to  fifth  order  was  used  and  appeared  to  yield  satisfactory  precision. 
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H(x)  =  H(x  )  +  ph  H'(x  ) 

°  0  2! 


H"(xo)  + 


+  <£*>_  H(V>(x  ). 
5!  0 


with 


H"(x  )  =  -2x  H'(x  ); 
o  o  o' 

H"'(xo)  =  (4x2q  -  2)  H'(xo),  etc. 


where 


ph  =  x  -  xQ  and  h  =  Xj-xo=x2-Xj  =  etc.  with  p  4^12. 

In  this  way,  values  for  H(x)  were  obtained  which  were  accurate  to  9 
or  10  decimal  places,  depending  upon  the  range  of  values  for  x  .  All  re¬ 
sults  were  stored  on  punchcards  for  use  in  the  main  program. 


Main  Program 

In  the  first  part  of  the  main  program,  four  different  types  of  isosceles 
triangles  were  considered,  characterized  by  : 


u 

■  \|  8/3, 

0~  no° , 

triangle  I;, 

u 

■  fn73. 

0~  146°, 

triangle  II; 

u 

*  2  . 

0  =  180°, 

triangle  III; 

u 

=  MT  . 

0  =  120°, 

triangle  IV. 

Each  triangle  was  considered  in  two  different  configurations  : 
i)  Rab  as  one  of  the  two  equal  sides  of  the  triangle; 
ii)  Rab  as  the  third  side  of  the  triangle. 

In  case  i)  the  results  for  <H'ab>  /0e2  and  <  H'^  ^/{Je2  were 
abbreviated  as  A  j,  Bj,  respectively,  whereas  in  case  ii)  these  results 
are  denoted  by  A2,  B2. 
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For  each  of  the  four  triangles,  and  for  each  value  of  the  parameter 


x  =  0Rab  the  following  quantities  were  calculated 

triangle  I 

V 

A-  ,  B 
l2  A1 

B  -  \  R-  • 

l2  1 

triangle  II 

AII1* 

AII  *  BII  * 
u2  ul 

BII2  ;  RII  • 

with 

R  = 
n 

2  A  ,  +  A  0 
nl  n2 

2  B  -  +  B  0 
nl  n2 

• 

The  ratio  Rn  gives  the  relative  threebody  component 


<Habc>  -  <H,abc> 
<H’abc> 


(o) 


for  a  triangle  of  type  n  . 


It  was  found  that  for  values  of  x  larger  than  3  the  results  become 
more  or  less  irregular.  This  is  due  to  the  fact  that  then  the  H-functions 
are  practically  equal  to  ^  2/n  and  the  H'  equal  to  zero;  partial  results 
may  become  as  small  as  one  in  the  ninth  or  tenth  decimal  place.  It  would, 
nonetheless,  have  been  possible  to  retain  high  accuracy  also  for  this  range 
of  values  for  x  by  using  a  double-precision  calculation  (to  22  decimal  pla¬ 
ces).  This  was,  however,  considered  unwise  since  it  would  have  prolon¬ 
ged  the  calculation  times  very  considerably. 


In  the  second  part  of  the  program,  two  more  types  of  triangles 
were  evaluated,  characterized  by  : 

u  =  1,  9  *  60°,  triangle  V 
u  *  2,  0  «  90°,  triangle  VI, 
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which  was  possible  without  principal  modifications  of  the  program.  The 
results  were  plotted  as  function  of  6,  and  proved  to  exhibit  regular  be¬ 
haviour. 

In  Fig.  1,  (section  B),  the  relative  threebody,  first-order  interac¬ 
tions  for  an  equilateral  triangle,  a  120° -symmetric  and  a  linear  symme¬ 
tric  array  are  plotted  as  a  function  of  x  =  |3R.  In  Table  III  detailed  results 
are  given  for  A  =  <H'ak  >  /p  e  for  the  two  different  configurations  of 

each  triangle  I,  II,  III,  IV,  V  and  VI  (Av  =  A  ),  as  well  as  the  twobody 

’  1  *9 

interaction  energy  By,  for  x  between  1.  5  and  3.d. 

Secondary  Programs 

The  foregoing  results,  transferred  to  punchcards,  were  applied  to 
a  determination  of  first-order,  threebody  interactions  between  triangles 
of  atoms  formed  by  a  central  atom  and  a  pair  of  nearest  neighbors  in  the 
face-centered  cubic  (fee)  and  hexagonal  close-packed  (hep)  structures. 

The  total  threebody  interaction  may  be  written  as 

N  (0**>  <EaAcp  and  N(Pe2)<EaAc- 

where  N  is  the  total  number  of  atoms  in  the  crystal,  £  is  the  electron 

charge,  and  E  is  the  total  threebody  interaction  for  a  central  atom 
a  A 

(a)  and  its  twelve  nearest  neighbors.  The  expressions  for  E  are  : 

<EaA»hcp  =  24  Av  +  18  <2  Arv1+Aiv2l  +  3  <2  Anij  *  Anij  > + 

*  12  (2  AyI  +  AyI^  )  +  3  (2  A,  +  A^  )  +  6  (2  Ajj  +  Aj^  )  ; 

the  roman  numerals  refer  to  the  different  types  of  triangles.  In  the  case 

of  an  equilateral  triangle  (type  V)  there  is,  of  course,  no  distinction  to 

be  made  between  Av  and  Av  ,  or  between  Bv  and  Bv  . 

V1  v2  V1  v2 
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<EaA>fcc  "  24  AV  +  18  (2  AIVl  + A!V2  >  +  3  (2  )  + 

+  12  (2  AyI^  +  AyI^  )  +  3  (2  Am^  +  A^)  +  6  (2  AJV^  +  A^)  . 

The  terms  have  been  arranged  in  such  a  way  that  the  first  four  are 
the  same  between  the  hep-  and  fcc-structures,  whereas  the  last  two  terms 
are  different. 

On  the  other  hand,  the  total  twobody  interaction  is 

N<Pe2)  Ea-b  *  with  Ea-b  =  6  BV* 

No  distinction  has  to  be  made  between  the  twobody  interactions  in 
the  hep-  and  fcc-lattices,  since,  as  far  as  a  central  atom  and  its  twelve 
nearest  neighbors  are  concerned,  these  interactions  are  the  same,  at 
the  same  value  of  x  =  p  R,  for  the  two  structures. 

The  following  quantities  were  calculated  as  functions  of  x  : 

^EaA^hcp  '  ^EaA  ^fec  ’  ^EaA^hcp  "  ^EaA^fcc  * 
together  with  ratios 

(E  ).  /  E  .  and  (E  ).  /  E  .  ; 

a  a  hep  '  a-b  aA  fee  '  a-b 

the  latter  quantities  give  the  relative  contributions  of  first-order,  three- 
body  interactions  for  the  two  lattices,  as  a  function  of  x  =  p  R.  Detailed 
results  are  given  in  Table  IV. 

A  final  subprogram  was  used  to  check  the  validity  of  the  Gaussian 
model  at  small  interatomic  distances.  To  this  end,  .the  repulsive  part  of 


the  modified  Buckingham-potential 
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£  (  6_  )  ea(i-R/Ro) 

1  -  (6/a)  a 

with  e  ,  £  and  Rq  given  parameters  for  neon,  argon,  krypton  and  xenon, 

was  compared  with  the  first-order  interaction  between  two  Gaussian  atoms, 
2 

(pe  )  By,  as  a  function  of  pR,  R  being  the  distance  between  the  two  atoms. 
Fig.  2  gives  the  results  of  this  comparison  (section  B). 
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IMPLICATIONS  OF  RESULTS  FOR  POSSIBLE  FUTURE  WORK 


The  evaluation  of  threebody  interactions  in  first-order  of  pertur¬ 
bation  theory  being  terminated,  the  next  step  consists  in  calculating 
the  threebody-second-order  (2-2)  effect.  This  term  is  on  one  side  mo¬ 
re  complicated  in  that  second-order  perturbation  theory  is  used,  but 
simpler  in  that  we  consider  here  the  interaction  between  a  diatomic 
unit  and  a  single  atom.  Electron  exchange  between  the  atom  and  the 
pair  is  consistently  neglected;  taking  it  into  account  implies  a  higher- 
order  correction  in  the  present  scheme. 

The  perturbation  hamiltonian  H'abc  takes  a  particularly  simple 
form  in  this  case.  Namely,  the  second -order  interaction  between  the 
two  atoms  which  form  the  diatomic  unit  (ab)  is  not  changed  by  the  pre¬ 
sence  of  the  third  atom  (c),  since  (c)  does  not  overlap  with  (ab).  Con¬ 
sequently,  for  the  evaluation  of  the  threebody  component,  we  may  res¬ 
trict  ourselves  to  the  part  H^abjc  of  H'abc  ,  giving  only  the  electro¬ 
static  interaction  between  the  pair  (ab)  and  the  third  atom  (c).  As  a  re¬ 
sult,  the  interaction  between  atoms  a_  and  b  does  not  occur,  which 
simplifies  the  calculations  considerably. 


Further,  the  summation  over  all  excited  states  of  the  three  atoms 
will  be  carried  out  by  replacing  the  resonance  denominators  by  an  ave¬ 
rage  excitation  energy,  which  leaves  only  matrix  elements  of  H'(ab)c 
and  (H'{ab)c)2  for  the  unperturbed  ground  state.  The  value  of  this  average 
excitation  energy  does  not  matter,  since,  it  drops  out  in  a  calculation  of 
the  relative  deviation  from  the  assumption  of  twobody  interactions.  The 
total  number  of  resulting  integrals  is  64,  mostly  of  three-center  type 
with  Gaussian  integrands,  except  for  factors  Vr^j  or  (^/r^jX  V1*^), 
where  i_,  j,  k,  _1_  denote  any  one  of  the  three  electrons  or  the  three 
nuclei,  with  i  /  j,  k  /  1.  These  integrals  are  in  part  direct,  in  part 
of  exchange  type;  the  exchange  involves  electrons  of  the  pair  (ab),  but 
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not  of  atom  c.  Further  stages  of  the  research  project,  such  as  the 
calculation  of  elastic  constants  (deviations  from  the  Cauchy-relations) 
will  then  be  devised  along  similar  lines,  pending  the  outcome  of  the 
present  computation. 


PUBLICATIONS  AND  CONFERENCES 


Three  Special  Technical  and  Scientific  Reports  were  issued  during 
the  first  year  of  the  Contract.  A  manuscript  entitled  "Systematic  Ana¬ 
lysis  of  Manybody  Interactions  in  Molecular  Crystals"  was  submitted  for 
publication  in  The  Physical  Review  (August  1961). 

Colloquium  talks  on  this  research  were  delivered  by  the  principal 
investigator  at  Ohio  State  University,  Columbus,  Ohio  (November,  1960), 
at  the  Max-Planck  Institut  fQr  Physik  und  Astrophysik,  Munich  (June, 

196  1)  and  at  the  National  Bureau  of  Standards,  Washington  D.  C.  (Octo¬ 
ber,  1961).  In  addition,  the  research  was  discussed  at  several  colloquia 
of  the  Battelle -Geneva  Institute. 


PERSONNEL  ON  THE  CONTRACT 


The  following  persons  were  employed  part-time  under  the  contract 
during  the  past  year  : 

1)  Dr  Laurens  Jansen,  principal  investigator; 

2)  Dr  Barbara  Kohin,  visiting  scientist,  the  Catholic  University 
of  America,  Washington  D.  C. ,  as  off  May  1961; 
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3)  Mr  Francis  Siegrist,  mathematics  student.  Federal  Institute 
of  Technology,  Zurich,  Switzerland,  during  August  1961; 

4)  Mr  Maurice  Marchais,  Mr.  Raymond  Lenoir,  Mrs  Jolan 
Targonski,  scientific  personnel  associated  with  the  digital 
computer  Bull  Gamma  AET. 


No  changes  in  research  policies  by  the  Contractor  have  occurred 
during  the  period  October  1,  1960  -  September  30,  1961. 


SUMMARY  OF  NUMBER  OF  MANHOURS  ON  CONTRACT 


The  total  numbers  of  manhours  expended  on  the  contract  for  the 
individuals  mentioned  in  the  previous  section  are  : 


-  Dr  Laurens  Jansen  :  1  000  hrs 

-  Dr  Barbara  Kohin  :  100  hrs 

-  Mr.  Francois  Siegrist  :  150  hrs 

-  Mr.  Raymond  Lenoir  :  125  hrs 

-  Mr.  Maurice  Marchais  :  250  hrs 

-  Mrs  Jolan  Targonski  :  300  hrs 


Not  listed  are  hours  expended  for  secretarial  services,  preparation 
and  processing  of  Technical  Reports  and  similar  items. 

No  important  property  was  acquired  during  the  contract  period  at 
direct  contract  expense. 


Geneva,  November  15,  1961 
LJ/gbs  -  Z-9219 


